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Abstract 

The 7/7 and 7/7 transition form factors are analyzed to leading- 
twist accuracy and next-to-leading order (NLO) of perturbative QCD. 
Using an rj — rf mixing scheme and all currently available experimental 
data the lowest Gegenbauer coefficients of the distribution amplitudes 
for the valence octet and singlet qq and the gluon-gluon Fock com- 
ponents are extracted. Predictions for the g*g*rf vertex function are 
presented. We also comment on the new BELLE results for the n^f 
transition form factor. 

1 Introduction 

The recent measurements of the photon to pseudoscalar meson transition 
form factors at large photon virtualities, Q 2 , by the BaBar collaboration 
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[1, 2] caused much excitement and renewed the interest in the theoretical 
description of these observables. Most surprising is the seemingly sharp rise 
of the 7T7 form factor with the photon virtuality, which is difficult to ac- 
commodate in fixed-order perturbative QCD. Power corrections to the usual 
leading-twist (collinear) approach [3] seem to be required. However, there 
is also a measurement of this form factor by the BELLE collaboration [4] 
which shed doubts on the BaBar data. At large Q 2 the results of the two 
measurements differ, the BELLE results are close to the theoretical expecta- 
tions from the leading-twist approach. Also the 777 and 7/7 form factor data 
behave as expected according to theoretical analyses [5, 6, 7] of the CLEO [8] 
and L3 data [9], measured at lower Q 2 than the BaBar data. Strong power 
corrections are not demanded in these analyses. Indeed a next-to-leading 
order (NLO) leading-twist analysis [6] is in reasonable agreement with the 
CLEO and L3 data. It is therefore tempting to reanalyze the 777 and 7/7 form 
factors to this accuracy along the lines presented in [6], taking into account 
the new BaBar data. It should be mentioned that the combined CLEO, L3 
and BaBar data on these two form factors have already been analyzed in 
other approaches. Thus, for instance, in [10] kj_ factorization is exploited 
or, in [11], light-front holographic QCD. In [12] the low Q 2 data have been 
studied within the non-local chiral quark model and in [13] the dispersive 
representation of the axial anomaly is used to derive an expression for the 
form factors that holds at all Q 2 . In [14] the anomaly sum rule has also been 
exploited for the analysis of the transition form factors. A combined analysis 
of the low and high Q 2 data has also been performed in [15]. 

The 77 and 77' mesons possess SU(3)f singlet and octet quark-antiquark 
Fock components and, additionally, two-gluon ones leaving aside higher Fock 
states. Leading-twist distribution amplitudes, 0, are associated with each of 
these Fock components. This feature leads, on the one hand, to the well- 
known flavor mixing and, on the other hand, as a further complication, to 
mixing of the qq singlet and the gluon-gluon (gg) distribution amplitudes un- 
der evolution. In the case of the transition form factors the two-gluon Fock 
components do not contribute to leading order (LO), they require higher 
orders of perturbative QCD. In our previous analysis [6] of the 777 and 7/7 
transition form factors to NLO leading-twist accuracy (see Fig. 1 for relevant 
Feynman graphs) the short lever arm provided by the Q 2 range of the CLEO 
and L3 data, the overall number of the data as well as the size of their errors, 
made it impossible to fix precisely even the lowest Gegenbauer coefficients of 
the distribution amplitudes. The resulting values were subject to substantial 
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errors. In particular we didn't find a clear signal for the gg distribution am- 
plitudes, a zero effect could not be excluded. Here, in this work we attempt 
a reanalysis of the form factors, using in addition to the CLEO and L3 data 
also the recent BaBar data [2]. With our analysis we want to demonstrate 
that, in contrast to the case of the pion, the r\ and 7/7 form factor data can 
easily be accommodated by a QCD calculation to NLO leading-twist accu- 
racy. A second goal of our analysis is the determination of the gg distribution 
amplitudes. On the strength of a larger set of data covering a substantially 
larger range of Q 2 , more precise results for the gluon distribution amplitudes 
than obtained in [6] can be expected and indeed, as it will turn out, is found. 
This is of utmost importance since the two-gluon Fock components play a 
role in many hard exclusive processes involving 77 and/or 7/ mesons which are 
nowadays accessible to experiment. Thus, for instance, the g*g*f]{f]') vertex 
contributes to decay processes such as T(1S) — > rj'X (see e.g. [16, 17]), to 
the hadronic production pp(p) — > rj'X and to meson pair production in the 
central region of proton-proton collisions [18]. The two-gluon Fock compo- 
nents may also matter in Xcj decays in pairs of rj or rj mesons and they may 
be attributed to the enhancement of some of the rj' channels in charmless B 
decays as compared to the corresponding pion channels (e.g. [19, 20, 21]). 
An example is set by the rj'K channels for which such an enhancement is 
experimentally observed [22, 23]. In this context the B — > 77(77') form factors 
are of importance which are also affected by the gg component of the 77 and 
rj mesons (e.g. [24]). 

The plan of the paper is the following: In Sect. 2 we recapitulate prop- 
erties of the qq and gg distribution amplitudes for the i] and 77' mesons, 
in particular their evolution behavior. The transition form factors to NLO 
leading-twist accuracy are presented in Sect. 3 while, in the next section, 
rj — rj mixing is briefly discussed. The CLEO [8] and BaBar [2] data on 
the transition form factors are analyzed for several scenarios in Sect. 5 and 
values for the second Gegenbauer coefficients of the three distribution ampli- 
tudes are extracted. We also shortly comment on the data for the time-like 
transition form factors obtained by BaBar [25]. Implications of the resulting 
distribution amplitudes for the g*g*P vertex are discussed in Sect. 6 and a 
brief comment on the ir / -f transition from factor is given in Sect. 7. As usual 
the paper ends with a summary. 



3 



7* 




Figure 1: Sample Feynman graphs contribution to the transition form factors 
to NLO. 

2 Properties of the distribution amplitudes 

When considering the evolution and mixing of quark and gluon states under 
evolution, it is convenient to choose SU(3) F octet and singlet combinations 
of quark-antiquark states [6] as valence Fock components of the pseudoscalar 
mesons, 

I QQs) = | {uu + dd — 2s s) / y/6) 

\qqx) = \(uu + dd + ss)/VS) (1) 

and the two-gluon state, \gg), which also possess flavor-singlet quantum num- 
bers and contributes to leading-twist order. Higher Fock components are 
neglected in our analysis since their contributions are power suppressed. It 
is to be stressed that the above states are partonic Fock components and 
not effective meson states or glueballs which are frequently considered in the 
treatment of 77 - rj' mixing. In collinear approximation and to leading-twist 
accuracy a distribution amplitude, <f>pi, is associated with each of the Fock 
components we consider 3 . They satisfy the symmetry relations (i = 1,8, 
P = 77,77') 

<f>Pi(x,fx F ) = Pi (l - x,/x F ) , <f> Pg (x,fx F ) = -<f>p g (l - x,/if) ■ (2) 

The quark distribution amplitudes are normalized to unity at any factoriza- 
tion scale ix F 

/ dx(f)pi(x,ix F ) = 1 (3) 
Jo 

3 A formal definition of the leading-twist distribution amplitudes in terms of particle- 
vacuum matrix elements of quark field operators or the gluon field strength tensor can be 
found in [6]. 
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while for the gluon distribution amplitude one has 

/ dx<f>p g (x,fx F ) = 0. (4) 
Jo 

Thus, there is no natural way to normalize the latter distribution amplitude, 
but its mixing with the quark distribution amplitude under evolution removes 
this ambiguity and, as shown in [6], the change in gluon normalization is 
reflected in the change of off-diagonal anomalous dimensions governing the 
evolution - for details, see below. 

The distribution amplitudes can be expanded upon the Gegenbauer poly- 
nomials 

<f>Pi{x,ti F ) = Qx(l-x)\l+ ]T a Pn (fi F )C,! /2 (2x-l) 



<f) Pg {x,fi F ) = X 2 {1 - Xf fl Pn(M C^\{2X - 1) , 

n=2,4,- 



(5) 



where only the terms for even n occur as a consequence of the symmetry 
relations (2). In terms of the expansion coefficients ap n the mixing of the 
quark singlet and gluon distribution amplitudes are expressed by 



4 n W = apnVo) - 



7i +) //9o 



+ pi ) a pl M 



a 9 Pn (ji F ) 



n u Pn 



' «s(/fo) 



7i + V/3o 



-+- a Pn 
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7i '/A 



,(6) 



where /iq is the initial scale of the evolution and (3o = 11/3N C — 2/3n/. The 
number of colors is denoted by N c and rif is the number of active flavors at 
the characteristic scale of the process. The coefficients of the eigenf unctions 
of the matrix evolution equation which introduces mixing between quark and 
gluon distribution amplitudes, are here denoted by Op^. The powers are 
the eigenvalues of the matrix of anomalous dimensions [26, 27] 
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where the LO elements of that matrix read (C F = (N 2 — 1)/(2N C )) 
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n > 2. 



(8) 



It is helpful for the understanding of the following to introduce the distinc- 
tion between the quantity Nf (= 3) which counts the (fixed valence) quark 
content of the meson's flavor- singlet combination and the above defined rif, 
the number of active flavors at some scale, which appears in the /3 functions 
and, as such, is connected to the running of the coupling constant. The 
parameters in (6) read 
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The evolution of the octet distribution amplitude is merely governed by 
and takes the simple form 



a 



{8) (n 



a, 



a, 



(10) 



To the order we are working, NLO evolution of the quark distribution 
amplitude should in principle be included. To this accuracy the Gegenbauer 
polynomials C T 3 / 2 are no longer eigenfunctions of the evolution kernel, i.e. 
their coefficients a % Pn do not evolve independently [28]. The impact of the 
NLO evolution on the transition form factors is expected to be small com- 
pared with the NLO corrections to the subprocess amplitudes [29]. Therefore 
we refrain from considering NLO evolution here. 

The projection of a collinear qq state onto a pseudoscalar state is achieved 
by replacing the quark and antiquark spinors by [30] 



T>i - C 

' afS,rs,kl ~ 




(11) 



where p denotes the meson momentum, while a (r, k) and j3 (s, I) represent 
Dirac (flavor, color) labels of the quark and antiquark, respectively. The 
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flavor matrices read 

C 8 = -^A 8 , Ci = -£=1/, (12) 

where \$ is one of the usual SU (3) p Gell-Mann matrices and 1/ is the 3x3 
unit matrix. The projector of a state of two (incoming) collinear gluons of 
color a and b and Lorentz labels \x and v carrying the respective momentum 
fractions x and 1 — x onto a pseudoscalar meson reads [6] 



where 



^ 2\jN f ,/n£-{ x(l - x) ' [6) 

ef" = ^ (14) 

(eoi23 = 1) describes the helicity-zero combination of the two transverse gluon 
polarization vectors. Taking for the meson momentum p = [p + ,0, 0j_], n in 
(14) is a light-like vector, n— [0,1, 0_J, which defines the plus component of a 
vector 4 , v + = n-v . The factor [x(l — x)]^ 1 appearing in the projector (13), is 
a consequence of the conversion of the potential of the gluon field occurring in 
a perturbative calculation of reactions involving two-gluon Fock components, 
into the gluon-field-strength operator whose vacuum-meson matrix elements 
define the gluon distribution amplitude, see [31]. The projectors are to be 
used in calculations of hard partonic amplitudes which are to be convoluted 
with either f l P / (2\/2N c )<j)p i or, for the gluons, fp/(2y/2N c )<f)p g . The decay 
constants, f P , are, as usual, defined by vacuum-meson matrix elements of 
flavor singlet or octet weak axial-vector currents 

(0 I 4s(0) I Pip)) =i/j*V (15) 

The singlet decay constants, fp, depend on the scale but the anomalous 
dimension controlling it is of order a s 2 [32]. In fact the evolution of f P1 is part 
of the NLO evolution of the singlet distribution amplitude, it represents the 
scale dependence of its first Gegenbauer coefficient. This is to be contrasted 
to the octet distribution amplitude for which the first Gegenbauer coefficient, 



4 Actually, in (14) instead of n any four vector can be used that has a non-zero minus 
and a vanishing transverse component. 
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do, is 1 at all scales. Thus, in harmony with the neglect of NLO evolution of 
the distribution amplitudes, the scale dependence of fpi is neglected too. 

Finally, we comment on different conventions of gluon distribution ampli- 
tude and singlet evolution to be found in literature. A detailed comparison 
is given in [6]. There are two important observations. First, as explained in 
detail in [6], a change of the definition of the gluon distribution amplitude 

<t>° Pg = a<t> Pg . (16) 

is (due to the fact that the physical quantity is independent of the choice 
of the convention) accompanied by a corresponding change in the hard scat- 
tering amplitude, i.e., in the projector of gg state, V 9 „ = 1/aV 9 , and 
(since the quark and gluon distribution amplitudes are connected through 
evolution) by a change of the off-diagonal anomalous dimensions 

7r* = -7„ w , 7^ = ^, (17) 

a 

or, equivalently, the evolution kernels V qg and V gq . The analogous relations 
hold for any change in the definition of the quark distribution amplitude, 
but since in contrast to the gluon distribution amplitude, its normalization 
is well defined, the unique definition of the quark distribution amplitude is 
encountered throughout the literature. 

Second, for historical reasons (see early work [27]) anomalous dimensions 
are in this field defined as in App. C of [6]: 

V AB (u,v) ® [weight B ](v)C B (2v - 1) = ^ B [weighty] (u)C A (2u - 1) , (18) 

where A,B<E {q, g} with corresponding notation C q = C^J 2 and C 9 = C^[l x 
while [weight J (u) = u(l — u) and [weighty] (u) = u 2 {l — u) 2 . On the other 
hand, the anomalous dimensions of the forward and non-forward parton dis- 
tributions (appearing in deeply inelastic scattering and deeply virtual Comp- 
ton scattering which are through operator product expansion closely con- 
nected to the meson transition form factors) are usually [33] defined as 5 

C A (2u - 1) <g> V AB (u,v) = - AB ln C B {2v - 1) . (19) 

5 Thc additional factor 1/2 present in [33] is compensated by factor 2 in their definition 
for the evolution kernel. 
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This definition leads to 




n{n + 3) 



36 




n{n + 3) 
36 



(20) 



Note that, as given in App. C of [6], 



n 



n(n + 3) 



(21) 



36 



Thus in this case the difference in off-diagonal anomalous dimensions is not 
due to different definition of gluon distribution amplitude but rather to dif- 
ferent defining equations of anomalous dimensions. When this, i.e., (20), is 
taken into account, one sees that our expressions are in agreement with [33], 
i.e., with the conventions used in deeply virtual processes. 

3 The transition form factors 

The 777 and 7/7 transition form factors can be represented as a sum of flavor- 
octet and singlet contributions 



where the singlet one also includes the gluon part. This decomposition is 
completely general. To NLO leading-twist accuracy the octet and singlet 
contributions can be written as [6] 



The meson and quark masses are neglected throughout. The a s correction 
to the quark parts, calculated in the MS scheme [34, 35], read 



F P7 (Q 2 ) = F P 8 7 (Q 2 ) + F P 1 7 (Q 2 ) 



(22) 
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while for the gluon contribution one has [6] 



1 1 — 3a; u 2 
K. g {x, Q 2 , hf) = —-In 2 a; In a; + In a; In — ^ . (25) 



The charge factors in (23) are 

n e u + e d ~ 2e 2 e 2 + e 2 , + e 2 
6 8 = , 6i = 7 = . (26) 

The charge of quark a in units of the positron charge is denoted by e a . The 
x -H- 1 — rr symmetry (2) of the distribution amplitudes has already been 
exploited in (23). As is well-known and can be seen from (5) - (8), the quark 
distribution amplitudes evolve into the asymptotic form 

0as = 6x(l x) (27) 

and the gluon one to zero for Q 2 — > oo. In this limit the transition form 
factors become _ 

2 f P + 2^2 ft 

3 (28) 

independently of the choice of the factorization scale. 

In our analysis we do not take into account power corrections as for 
instance may arise from higher-twist effects 6 , from meson masses or from 
quark transverse degrees of freedom. The latter seem to be rather small for 
the form factors of interest, in fact much smaller than for the 7T7 transition 
form factor as is shown for instance in [10]. For further comments concerning 
power corrections see Sect. 5. 




4 r] — rj mixing 

In order to reduce the number of independent distribution amplitudes we 
follow [37] and assume meson independence of the distribution amplitudes, 
i.e. 

4>Pi = 4>i, 4>Pg = 4>g (29) 

6 It can be shown that there is no twist-3 correction to the transition form factor. 
Possible twist-4 and twist-6 corrections have been discussed for the case of the form 
factor recently [36]. 
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as we did in [6]. Hence, the mixing behavior of the valence Fock components 
and the particle dependence solely resides in the decay constants. Since in 
hard processes only small spatial quark-antiquark separations are probed, 
this assumption is sufficiently plausible - the decay constants play the role of 
wave functions at the origin of configuration space. 

For the decay constants we use the general parametrization [32, 37] 

f'l = h c °s #s , f'l = ~fi sin 0i , 

4 = / 8 sin0 8 , fl f = fi cos 9± . (30) 

In [37, 38] it has been observed that rj — rj' mixing is particularly simple in 
the quark-flavor basis defined by the transform 

| rj q ) = cos y? | if + sin y? | rj') , 

| rj s ) = - sin y? | if + cos y? | rj') , (31) 

of the physical meson states. In this basis which is supposed to separate 
strange and non-strange contributions, the mixing behavior of the decay 
constants is controlled by the angles ip q and ip s , defined analogously to (30). 
It turned out from phenomenology that these angles practically fall together 
with the state mixing angle <p, i.e. that we can write 

/* = f q cosip, /*=-/ s siny?, 

= f Q smtp, f°,= / s cosy?. (32) 

This observation is supported by a QCD sum rule study [39]. A recent 
lattice QCD study [40] is also not in conflict with it. The occurrence of only 
one mixing angle in this basis is a consequence of the smallness of OZI rule 
violations which amount to only a few percent and can safely be neglected 
in most cases. SU (3) f symmetry, on the other hand, is broken at the level 
of 10 — 20%. The assumption (p q = y? = (p s is essential for the quark-flavor 
mixing scheme. 

Working in the quark-flavor basis and exploiting the divergences of the 
axial-vector currents - which embody the axial-vector anomaly - the mixing 
parameters in the quark-flavor mixing scheme can be expressed in terms of 
the masses of the physical mesons [41], e.g. 

2(M* -M^M^-M^) ■ [ ] 



sin if 



\ 
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Ref. 


hi U 


#8 




0i 




1 


[41] 


1.19 


-19.4° 


1.10 


-6.8° 


41.4° 


2 


[37, 38] 


1.26 


-21.2° 


1.17 


-9.2° 


39.3° 


3 


[42] 


1.51 


-23.8° 


1.29 


-2.4° 


40.7° 



Table 1: Decay constants in the singlet-octet basis and the mixing angle in 
the quark-flavor basis. The value of tp quoted in the last line for is the average 
of ip g and tp s determined in [42] . 

Evaluation of the mixing angle provides tp = 41.4°. Using for the decay 
constant f q in the quark-flavor basis the SU(3)p symmetry result f q = f n , 
one finds for the strange decay constant in that basis 

2(M2 -2M| + M2 )(2M|. -M2 -M2)- 1 > 

Transforming these results to the octet-singlet basis, one obtains the results 
for the mixing parameters (30) that are quoted in Tab. I. We repeat it is 
convenient to work in the latter basis because in it evolution which is an 
important ingredient to our analysis, is simple. In [37, 38] the mixing pa- 
rameters have been determined phenomenologically using tp q = tp = tp s but 
allowing for higher-orders flavor symmetry breaking effects. The results, ob- 
tained from the analysis of a number of processes involving 77 and 7/ mesons, 
are also quoted in Tab. 1. For a discussion of uncertainties see [37]. There 
are a few more analyses, e.g. [39, 42, 43, 44], in which the decay constants 
(30) have been determined. In general the results are close to those ob- 
tained in [37, 38], for a comparison see [41]. The largest deviations from 
the mixing parameters given in [37, 38] has been reported in [42], see Tab. I. 
This phenomenological analysis has been performed along the lines described 
in [37, 38] considering however only a subset of the processes investigated 
therein but, if at disposal, exploiting more recent data. In other papers only 
the mixing angle tp has been determined, e.g. [45, 46]. Within occasionally 
large errors the values for it (39 — 42°) agree with the ones found in [37, 38]. 

Frequently, mixing of the i] and rf is studied starting from three basis 
states, r] 8 ,rji ( or r) q ,r) s ) and a gluonic state r) g . Despite this the values for the 
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angle ip controlling rj-rj' mixing obtained in such analyses agree reasonably 
well with the above quoted ones. However, these mixing schemes rely on 
the existence of a rather light pseudoscalar glueball for which there is no 
evidence, see the review [47] . In any case, it would be a misinterpretation to 
identify the gluon-gluon state we consider here with the r) g . Our gluon-gluon 
state is a partonic Fock component of the rj and r] f . 



5 Analysis of the form factor data 

Let us now turn to the extraction of the various distribution amplitudes or 
rather their Gegenbauer coefficients from the data on the transition form fac- 
tors [2, 8, 9]. Inserting the Gegenbauer expansion of the distribution ampli- 
tudes (5) into (23) one can perform the integration term by term analytically. 
The result is 



Q 2 F, 



P7 



= a%((i F ) 1- 
+ af 2 (fi F ) 

20 a a (fi R 



Ti 



5 a s (fi R ) /59 _ 5 Q\ 

3 7T V 7 2 6 n fj? F ) 

5 ^(^ / 10487 91 

3 7T V 4500 75 11 
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1296 
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108"" <i 2 -> 



10125 



7 , Q\ 
675 ln W 



H F < 
+ ... 



where we introduced the abbreviation (n = 0, 2, 4, . . ., Oq = 1) 



(35) 



l °Pny 



M = \ -\fp<^ F ) + 2V2f P a l M 



(36) 



Before we analyze the data a comment is in order: Inspection of (35) 
reveals the familiar result that only due the admittedly mild logarithmic Q 2 
dependence generated by the evolution and the NLO corrections one can in 
principle discriminate among the Gegenbauer coefficients of different orders. 
However, even with the large range of Q 2 in which form factor data are 
available now, this logarithmic Q 2 dependence is in practice insufficient to 
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KJ.KJO 


/4 = QV2 


36.9 


-0.01 ±0.02 


-0.08 ±0.01 


10 ±4 


0.03 




37.5 


-0.01 ±0.02 


-0.07 ±0.01 


6±3 


0.03 


mixing [42] 


45.0 


0.05 ±0.02 


-0.16 ±0.01 


11 ± 5 


0.10 


4 = a§ 


49.8 


-0.11 ±0.01 


-0.11 ±0.01 


21 ±5 


0.0 


a2 = Og, mix. [42] 


240 


-0.16 ±0.01 


-0.16 ±0.01 


19 ±4 


0.0 



Table 2: Gegenbauer coefficients at /io = 1 GeV fitted to the data from [2] 
and [8] for Q 2 > 2 GeV 2 (22 and 18 data points, respectively). Except stated 
otherwise, the standard setting described in the text (with \i F = /i R = Q), 
and the mixing parameters of [37] are used (see Tab. 1). Eq. (41) is probed 
at no — 1 GeV. 

allow for an extraction of more than one Gegenbauer coefficient for each of 
the distribution amplitudes [48]. In view of this problem, we are forced to 
truncate the Gegenbauer series at n = 2. To leading-twist accuracy the 
higher Gegenbauer coefficients are not suppressed as is the case for other 
approaches in which power corrections, accumulated in the soft end-point 
regions x — > or 1, are taken into account [10, 36]. Therefore the n = 2 
coefficients we are going to determine below, suffer from a truncation error; 
they are to be viewed as effective parameters which are contaminated by 
higher order Gegenbauer coefficients. 

5.1 Fits 

Except stated otherwise we employ the following specifications in our fits: As 
the minimum value of Q 2 used in the fits we take 2 GeV 2 and for the initial 
scale of the evolution we choose /x = 1 GeV. For a s we use the two-loop 
expression with four flavors (rif = 4) and = 319 MeV [49]. For the 
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Figure 2: The 777 (left) and 7/7 (right) transition form factors scaled by Q 2 . 
Data taken from [2, 8, 9]. The solid lines represents our default fit (see Tab. 
2). The dotted lines are the LO asymptotic results (28). 

factorization and renormalization scales we adopt the frequently used choice 
fJ'F — I^r — Q which conveniently avoids the ln(Q 2 /fip) terms in (35) and, 
hence, a contingent resummation of these logarithms [51, 50]. In Sect. 5.1.1 
we will comment on other choices for these scales. 

Using this standard setting together with the mixing parameters derived 
in [37], we fit the Gegenbauer coefficients of order 2 to the data 7 8 9 [2, 8]. 
The result of this fit, termed the default fit in the following, is quoted in 
Tab. 2 and shown in Fig. 2. For comparison the results of fits to only the 
CLEO data [8] and only the BaBar data [2] are also shown in Tab. 2. The 
latter fit can be regarded as a change of the minimal value of Q 2 for which 
data are taken into account in the fits. One notices that the three sets of 
parameters for the quark distribution amplitudes agree quite well with each 
other. Deviations of a little more than la are seen for af. The effect of the 
Babar data in combination with the CLEO one results in a reduction of the 
parameter errors and a more precise determination of a%. In contrast to our 
previous work where we have had at disposal only the CLEO data [8] and 
have chosen fi R = Q/\f2 and fxp = Q, the gluon distribution amplitude is 

7 The signs of the transition form factors are not measured. 

8 The CLEO collaboration [8] measured the form factors for various rj(r]') decay chan- 
nels. We take into account all of them. In cases where for a given value of Q 2 there are 
several values of the form factor we use their error-weighted average. 

9 Above Q 2 — 2 GeV there is only one 7/7 data point from [9] with a large error. It 
has no bearing on our fits and will therefore not mentioned explicitly in the following. 
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Figure 3: The distribution amplitudes specified by the default fit at the initial 
scale. 

not compatible with zero now. The reasonable agreement of the three fits 
demonstrates the consistency of the CLEO [8] and BaBar data [2]. The x 2 of 
the default fit is very good given that all together 40 data points are included 
in the fit. In Fig. 2 the fit is compared to experiment. 

The distribution amplitudes corresponding to the default fit (see Tab. 2) 
are shown in Fig. 3. Those for the octet and singlet qq components are close to 
the asymptotic form of a meson distribution amplitude. They are symmetric 
around x = 1/2 while the gluon distribution amplitude is antisymmetric. 
We repeat that the distribution amplitudes are to be considered as effective 
ones since the parameters ci2 are contaminated by higher order Gegenbauer 
coefficients. 

5.1.1 Scale dependence 

Let us discuss the choices of the scales in more detail. In principle, renormal- 
ization and factorization scales can be chosen independently. In this work 
we have chosen fx R = up = Q for the default fit but we also investigate other 
choices of the scales in order to learn about the theoretical uncertainties of 
our results. 

First, we discuss our factorization scale choice. As elaborated in [50], the 
dependence on fip can be cancelled order by order in a s by performing the re- 
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summation of (a s In fJ. 2 F /Q 2 ) n terms 10 , i.e. by performing the evolution of the 
hard scattering part from \x F to Q the same way the logs (a s In fi^/ fi 2 F ) n are 
resummed in the distribution amplitude by employing the evolution equation. 
Occasionally, however, one evolves the distribution amplitude(s) from /i to 
/if, and does not resum the logs in the hard-scattering amplitude, thus caus- 
ing a residual dependence on \i F even after the convolution with distribution 
amplitude(s). But the choice fi F = Q corresponds to evolving distribution 
amplitude(s) all the way from fio to Q and removing all In fip/Q 2 terms from 
the hard-scattering part. This result is equivalent to evolution of distribution 
amplitude from /i to any ji F and the resummation of (a s \n n 2 F / Q 2 ) n terms 
in the hard-scattering part and in this sense the choice fi F = Q appears 
convenient. 

For the renormalization scale the usual choice is fi^ = Q. On the basis of 
a next-next-to-leading order calculation of the pion transition form factor it 
has been argued in [51, 50] that another possible choice of jiR is the square 
root of the average virtuality instead of the characteristic scale of the process, 
Q. For the transition form factors the average virtuality is Q 2 /2. Hence, 
in order to explore the renormalization scale dependence of the resulting 
distribution amplitudes we also utilize this choice. As an inspection of Tab. 
2 reveals the dependence of the fit on the renormalization scale is rather 
strong; the change of (j,r from Q to Q/y/2 results in substantial change 
of the parameters although both the fits are of similar quality. Thus, the 
theoretical uncertainties of our results are larger than indicated by the errors 
of the fitted Gegenbauer coefficients. 

It is important to note that our NLO calculation is in fact the LO calcula- 
tion in a s and it is a well known fact that in order to stabilize the dependence 
on ixr the NLO QCD corrections i.e., NNLO corrections to the transition form 
factor, should be included. At NNLO the presence of a 2 In \j? R j Q 2 terms sta- 
bilizes the dependence on /i R and all predictions fall relatively close (see, for 
example, [52] for discussion on that point). Concerning fits, one would thus 
expect that the inclusion of NNLO would considerably decrease the variation 
of the obtained Gegenbauer coefficients with /xr. Without this stabilizing ef- 
fect of NNLO we are left with the variation illustrated in Tab. 2. 

Finally, we test the strength of the residual dependence on the factoriza- 

10 The resummation/evolution in the hard-scattering part should be performed from 
(ip to the characteristic scale of the process which appears in the logs - in our case Q. 
One might say that the question then only remains what is the characteristic scale of the 
particular process but this is more of a problem for the two scale processes. 
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Scales 


1 9 

a 2 — a 2 


a% - a\ 


al — a 9 2 




0.371 


0.248 


0.057 


/4 = Q 2 /2 


0.607 


0.228 


0.057 


A = /4 = 372 


0.760 


0.202 


0.058 



Table 3: Correlation coefficients obtained in the fits to the CLEO [8] and 
BaBar [2] data for various choices of the scales. Except of the scales the 
standard setting and the mixing parameters of [37] are used. 

tion if the resummation of the (a s In n^/ fi 2 F ) n terms is ignored. In order to 
do so we perform a fit with the choice Hr = fj,p = Q/y/2. From the results, 
quoted in Tab. 2, one sees that the fit mildly depends on /ip. The Gegen- 
bauer coefficients al and al hardly change. For a 2 the scale dependence is 
a bit more pronounced although the values of a 2 agree within errors for the 
two fits with different [ip but the same renormalization scale. 

5.1.2 Parameter correlations and evolution 

Due to their mixing under evolution we expect the strongest correlation be- 
tween the Gegenbauer coefficients al and a 2 . Inspection of Fig. 4 reveals 
that the correlation is particularly strong for the fit to just the CLEO data 
but becomes milder for the fit to the combined CLEO and BaBar data. This 
goes parallel with a reduction of the parameter errors. The strength of the 
correlation between al and a 2 depends on the chosen factorization and renor- 
malization scales; it is smallest for the standard setting, see Tab. 3 where we 
compile the correlation coefficients for the three fits to the CLEO and BaBar 
data. On the other hand, the correlation between al and al is milder while 
a 2 and al are nearly uncorrelated. The origin of the correlation between a| 
and the other Gegenbauer coefficients lies in r\ — rj' mixing. The form factor 
F m is dominated by the octet contribution (af) while F^^ is mainly fed by 
the singlet one (al,a 2 ). 

The evolution of the Gegenbauer coefficients al and a 2 is shown on the 
right hand side of Fig. 4. The coefficients decrease relatively fast from the 
initial scale, up to about 10 GeV 2 and after that the approach to zero is 
rather slow. This behavior is a consequence of the properties of the logarithm 
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H 2 F [GeV 2 ] 



Figure 4: Correlation between a\ and af (left) and their evolution (right) for 
the standard setting. The la x 2 -contours are shown for fits to the data from 
only CLEO [8], only BaBar [2] and to both data sets. The shaded bands 
indicate the errors of the coefficients. 



contained in a s (fi F ) (see (6)): for fi F 3> //q (in fact for ji F > 3 GeV) the 
derivative of the logarithm is very small. The Gegenbauer coefficients of 
the octet distribution amplitude evolve similarly. This flat behavior of the 
logarithm for large scales is also partially responsible for the fact that even 
in a NLO calculation one can determine only one Gegenbauer coefficient for 
each distribution amplitude or, more precisely, one linear combination of 
them. 



5.2 OZI rule constraints 

The valence Fock components of the SU(3)i? octet-singlet basis can be trans- 
formed to those of the quark-flavor basis [6] 



\Vq) - 

\Vs) 



2y/2N c 
fs 



2y/2N c 



<f) q (x,ii F ) \qq) + 4>op P (x,fi F ) \ss) + J- <j> g (x,n F ) \gg) 



<Po PP {x, Hf) \qq ) + <j>s{x, n F ) \ ss ) + -j=(pg(x, fi F ) \gg) 



V3 



(37) 



where qq is short for the combination (uu + dd) / \/2 and 



i 9 = i(0 8 + 20O, ^ = ^(20 8 + 0i), O pp = ^(0i-0s). (38) 
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The new decay constants are related to / 8 and fi by 



/, = pR - fa , /. = pfi ~ fi ■ (39) 

We see that in (37) the ss (qq) Fock component appears in the r] q (r] s ) states. 
These respective opposite Fock components lead to violations of the OZI 
rule if they were not suppressed. Hence, the consistency of the quark-flavor 
mixing scheme and the perturbative approach to hard exclusive processes 
requires that 

l%^liff)J « 1 (40) 

holds for any values of x at least for a limited range of the factorization scale. 
Asymptotically where all distribution amplitudes evolve in 0as? 0o PP is zero 
anyway. To a sufficient degree of accuracy (40) can be replaced by 

^l«2-«2l«l (41) 

Indeed the default fit meets (41) at the initial scale and, as can readily 
be checked, at all larger factorization scales: ^ \a\ — a|| <0.03. Hence, no 
substantial violations of the OZI rule follow from the distribution amplitudes 
specified by the default fit and shown in Figs. 3. Also the fits using different 
choices of fx R and fx F as well as those to just the CLEO or BaBar data satisfy 

(41). ( ^ 

Using the 77 — 77' mixing parameters determined in [42] (see Tab. 1) which 
markedly differ from those given in [37, 38], one also arrives at a reasonable fit 
to the form factor data with regard to x 2 - However, in this case a% is positive 
leading to rather large OZI rule violations. The mixing parameters given in 
[42] are therefore to be employed with reservation. It is to be mentioned that 
the work [42] has also been criticized in [13] on the basis of an analysis of 
the transition form factors with a [/^(l^-aiiomaly sum rule. The theoretical 
set of mixing parameters discussed in [41] is intermediate between [37, 38] 
and [42]. The quality of the fit to the form factor data is similar to other 
fits but the difference I rather large although smaller than for 

the fit using the mixing parameters given in [42]. Thus, with regard to the 
strength of OZI rule violation the mixing parameters given in [37, 38] seem 
to be favored. Ideally one should fit the mixing parameters together with 
the lowest Gegenbauer coefficients to the data. Unfortunately such a multi- 
parameter fit does not lead to a reasonable solution, there are extremely 
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strong correlations among the parameters, often large violations of the OZI 
rule and a covariance matrix that is not positive definite. Thus, we refrain 
from discussing such fits. 

In order to avoid large violations of the OZI rule one may follow a sugges- 
tion made in [24] and assume a\ = a% at the initial scale. Although evolution 
generates some violations of the OZI rule with increasing scale they are always 
sufficiently small. The fit assuming a\ = a\ and using the mixing parameters 
of [37, 38] is still of reasonable quality. Although x 2 is somewhat larger, the 
fit parameters are similar to those obtained from the three-parameter fit, see 
Tab. 2. The difficulties with the OZI rule of the mixing parameters deter- 
mined in [42] is corroborated by the analogous fit with a\(n Q ) = a%(/j, ). This 
fit fails badly, X 2 is 240. 

5.3 Comparison to other results 

Ali and Parkhomenko [16] have performed an analysis of the rj' energy spec- 
trum in the inclusive decay T(1S) — > ggg* — > rj'X [53] in order to constrain 
the r/'-meson distribution amplitude. At an intermediate step of the analysis 
of the Y(IS') — > rj'X energy spectrum the g* — > vj ' g transition form factor 
is to be calculated (see Sect. 6). From the high end of the rj' meson energy 
spectrum Ali and Parkhomenko determine the Gegenbauer coefficients a\ 
and af which are in agreement with our results within very large errors. In 
order to calculate the energy spectrum also at low and even negative gluon 
virtualities Ali and Parkhomenko introduce a positivity constraint for the v( g 
transition form factor which is achieved by choosing \j? F = [i 2 R = \Q 2 \ + m 2 ,. 
This constraint significantly reduce the allowed range for the Gegenbauer 
coefficients and, in combination with our previous result [6] they obtain 

4 = -0.08 ±0.03, af = 6.5±2.7. (42) 

at the initial scale Hq = 1 GeV. While the value for a\ agrees with our default 
result within errors, our result for the coefficient af is larger as a consequence 
of the BaBar data which were not available to the authors of [16]. The role 
of the positivity constraint used in [16] remains to be understood. 

The lowest Gegenbauer coefficient of the octet distribution amplitude has 
been calculated with the help of QCD sum rules [54]. A value of about 0.2 
with a large uncertainty has been obtained for the coefficient a\. In this cal- 
culation the distribution amplitude is normalized to f% = fn- For the larger 
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octet decay constant (see Tab. 1) we are using the value of a\ is expected 
to be somewhat smaller, say, about 0.16. Given the large uncertainty of the 
sum rule result this value is not in conflict with our default fit. However this 
large face value of the sum rule result can only be reconciled with the form 
factor data if one allows for an additional power correction (or for higher 
Gegenbauer coefficients) which compensate the positive a% to a large extent. 
In order to examine this possibility we extract the octet contribution to the 
transition form factors using (22), (23) and (30) 

cos e x F m + sinflii^ 
mi ~ cos(9 8 -9 1 ) [ } 

and perform a fit analogously to the ones described above but allowing for an 
additional power correction c 8 /Q 4 to it. Keeping a\ = 0.16 fixed and using 
the mixing parameters determined in [37, 38], we obtain c§ = —0.08 ± 0.01 
from the fit to the data on F mi . For the 15 data points the minimum x 2 
is 21 which is somewhat larger than our best result (x 2 = 15.3 for a\ = 
—0.01 ± 0.02, c 8 = 0) but still tolerable. Freeing also a very good fit is 
obtained (a 8 2 = 0.06 ± 0.05, c 8 = -0.04 ± 0.02, X 2 = 16.3). The latter fit is 
shown in Fig. 5 and compared to experiment. In this figure also the results 
of the default fit for the octet and singlet form factors 

_ cos flsify-y - sin 9 8 F m 

*™~ COS (0 8 -0O 1 ] 

are displayed and compared with experiment. 



5.4 A comment on the time-like data 

The BaBar collaboration [25] has measured the 777 and 7/7 transition form 
factors at s = 112 GeV 2 . We do not include these data in our analysis since 
the theoretical treatment of form factors in the time-like region is non-trivial 
and not well understood. In the time-like region the form factors are no 
longer real and there are subtleties regarding the analytic continuation from 
the space- to the time- like region [55] (and references therein). 

Compared with the naive expectation that at s = 112 GeV 2 the form 
factor should be close to the asymptotic prediction the Babar result for the 777 
form factor is about 2a too large. This discrepancy is likely not a consequence 
of the description of 77 — rf mixing. The quark-flavor mixing scheme is on 
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Figure 5: The octet and singlet components of the transition form factors 
evaluated with the mixing parameters given in [37, 38]. Data taken from 
[2, 8]. The solid and dotted lines represent our default fit (see Tab. 2) and 
the asymptotic results, respectively. The dashed line is the fit with a power 
correction (a| = 0.06 ± 0.05, c 8 = -0.04 ± 0.02). 

sound theoretical grounds and is phenomenologically well established in a 
large variety of processes. A 2a discrepancy for a single data point cannot 
discard this mixing scheme. In contrast to the case of the r] the 7/7 form 
factor at s — 112 GeV 2 measured by the BaBar collaboration [25], is in 
agreement with the naive expectation. 

6 The g*g*P vertex 

We are now in the position to repeat the evaluation of the g*g*P form factors 
we performed in [6]. With the results on the gluon distribution amplitude 
obtained from the present analysis of the P7 transition form factors we be- 
lieve to have more precise results for the g*g*P form factors now. On the 
importance of these form factors we have already commented in the intro- 
duction. 

The gluonic vertex is written analogously to the electromagnetic one as 
T% = % F Pg *(Q 2 , to) 5 ab q la q 2 p , (45) 

where q% and q 2 now denote the momenta of the gluons and a and b label the 
color of the gluon. We consider space-like gluon virtualities for simplicity. It 
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Figure 6: Sample lowest order Feynman graphs for the g*g* — > qq and g*g* — > 
gg subprocesses to the g*g* — > P transition form factors. 

is convenient to introduce an average virtuality and an asymmetry parameter: 

1 



Q 2 = — 2 {ql + ql), 



id 



Qi + <?1 



(46) 



The values of u> range from —1 to 1, but, due to Bose symmetry, the transition 
form factor is symmetric in this variable: Fp g *(Q 2 ,u>) = Fp g *(Q 2 , —oj). 

According to [6] the Pg* transition form factor to leading-twist accuracy 
and lowest order of a s is to be calculated from Feynman graphs of which 
examples are shown in Fig. 6. It reads 



f 1 

F Pg *(Q 2 ,u) = 4vrct s (/V ' '" 



Q 2 N c 



+ 0(a 2 s ), (47) 



where 



i 2 i 

dx 4>i(x, /i F ) - — . , 
o 1 — ur(l — 2x) z 



l-2x 



XX 



l-u 2 (l-2x) 



(48) 



There are a few more LO Feynman graphs involving triple or quadruple gluon 
vertices. However, they do not contribute when contracted with the qq or gg 
projectors, (11) and (13). The integrals in (48) can be worked out by inserting 
the Gegenbauer expansions (5) into it and performing the integrals term by 
term analytically analogously to (35). The results are given in [6] and are not 
repeated here. Using the Gegenbauer coefficients of the default fit, we can 
readily evaluate the g*g* — > v( transition form factor. The results, including 
the la error band, are shown in Fig. 7 for two values of Q 2 = fi 2 R = n 2 F . As 
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Figure 7: Predictions for the g*g*rj' form factor evaluated from the default 
solution with fi 2 R = fi 2 F = Q 2 . The shaded bands indicate the la uncertainty 
of the predictions. 

compared to our previous results [6] the error bands are markedly narrower 
while the central values do not differ much. 

The ratio of the rjg* and rfg* form factors is given by 

Fng./Fw = fl/fl, = -tan0!. (49) 

As for the 777 and r/7 transition form factors and in order to be consistent 
with that analysis power corrections are neglected here as well. For the 
g*g* r l{ r l') vertex function, in particular for u — > 0, any power corrections as 
for instance quark transverse momenta or meson mass corrections, are small 
since the vertex function is not end-point sensitive in this kinematic limit. 
Mass corrections to the g*g*rj' vertex function have been estimated by Ali 
and Parkhomenko [16, 17]. 

7 The 777 transition form factor 

For comparison we show in Fig. 8 the data on the 717 transition form factor 
[2, 4, 8]. It is important to realize the dramatic difference between the BaBar 
data [2] on the tcj form factor and the T](r]')j data. In contrast to the latter 
a leading-twist analysis of the Try form factor to fixed order of perturbative 
QCD fails because the data do not seem to respect dimensional scaling. 
Power corrections seem to be demanded by the data as well as a positive 
value of the 2nd Gegenbauer coefficient of the pion distribution amplitude, 
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Figure 8: The 7T7 transition from factor. Data taken from [2, 4, 8]. The 
solid (dashed, dotted) line represents a fit with a^O^o) = —0.02 ([29], LO 
asymptotic result). 

see for instance [10, 14, 36]. This difference between the and the other 
form factors implies a strong breaking of flavor symmetry which has never 
been observed in the sector of pseudoscalar mesons before. 

On the other hand, the BELLE data [4] do not show this sharp rise with 
Q 2 , they lie systematically below the BaBar data. They do not exceed the 
asymptotic result for the ttj form factor and are close to the behavior of 
the NLO leading-twist approach as is the case of the other transition form 
factors. An example of such a fit to the data from [8] and [4] is shown in Fig. 
8 (with a£Oo) = -0.02 ± 0.02 and x 2 = 34.9 for 28 data points). Although 
the x 2 °f this fit is reasonable the Q 2 dependence of the fit seems to be too 
flat as compared to the BELLE data which may be viewed as a hint at lacking 
power corrections. For instance, a result obtained with kj_ factorization which 
encodes corrections of order (k\/Q 2 ), is in good agreement with the CLEO 
and BELLE data [29]. We finally emphasize that the BELLE data imply 
only mild violations of flavor symmetry. A more detailed comparison of 
the BaBar and BELLE data with theoretical models can be found in [60]. 
Although the BELLE data seem to be favored against the BaBar data with 
regard to the standard theoretical concepts, an understanding of the origin 
of the discrepancy between the two measurements is required. 
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8 Summary 



We have analyzed the 777 and 7/7 transition form factors within a collinear 
factorization approach to leading-twist accuracy and NLO of perturbative 
QCD. The analysis of the data [2, 8] allowed for an extraction of the lowest 
Gegenbauer coefficients of the qq flavor-octet and singlet distribution ampli- 
tudes as well as that of the glue-glue distribution amplitude. The Gegenbauer 
coefficients are better determined now and have smaller errors than those ex- 
tracted in [6]. Our default result for /j — 1 GeV is (see Tab. 2): 

al = -0.05 ±0.02, a\ = -0.12 ±0.01, af = 19 ± 5 . (50) 

There are a number of uncertainties of this result. First there is the uncer- 
tainty due to the chosen renormalization scale. Next, despite the large range 
of Q 2 covered by data now one still cannot determine more than the lowest 
Gegenbauer coefficients of the three distribution amplitudes. The coefficients 
we quote are to be regarded as effective parameters which are contaminated 
by higher order Gegenbauer coefficients. In fact, according to the discussion 
at the end of Sect. 5.1.2, it seems impossible to extract more information 
from the transition form factors than the lowest Gegenbauer coefficients. In 
order to determine more Gegenbauer coefficients additional processes have 
to be analyzed. Finally, given the quality of the present data, it is still not 
possible to discriminate between the logarithmic Q 2 dependence generated 
by the evolution and the NLO corrections, and power corrections (see the 
discussion in Sect. 5.3). We have neglected power correction in our analy- 
sis. Allowing for power corrections in the form factor analysis one may for 
instance find positive qq Gegenbauer coefficients as is the case for the pion 
distribution amplitude (e.g. [36]). 

As an application of our results for the Gegenbauer coefficients we calcu- 
lated the g*g*f]' vertex function. Another processes in which the gluon-gluon 
Fock components of the 77 and rj' mesons may play an important role, are the 
Xcj (J = 0,2) decays into pairs of rj and/or rj' mesons. Here, the c-quark 
mass is considered to be large enough to allow for a perturbative treatment 
of these decays. An explicit calculation of the perturbative contribution to 
the Xcj decays taking into account the two-gluon Fock components, is te- 
dious, many Feynman graphs contribute even to lowest order [56] . Moreover, 
there is another complication. As is shown in [57] the next higher Fock state, 
ccg, of the Xcj, the so-called color-octet contribution [58], is also to be taken 
into account since it scales with the same power of the hard scale, m c , as 
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the cc contribution. With regard to these complications whose detailed cal- 
culation is very time-consuming, we will not attempt a complete analysis 
of these decay processes; this is beyond the scope of the present paper. A 
statement whether or not our gg distribution amplitudes are in conflict with 
the peculiar features of the Xcj decays [46, 49, 59] is therefore premature. 

This work was supported in part by the BMBF under the contract No. 
06RY9191 and in part by Croatian Ministry of Science, Education and Sport 
under the contract no. 098-0982930-2864. 
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